Performance of prototypes for the ALICE electromagnetic calorimeter 
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The performance of prototypes for the ALICE electromagnetic sampling calorimeter has been 
studied in test beam measurements at FNA L and CE RN. A 4 x 4 array of final design modules 
showed an energy resolution of about 11% / E{GeV) © 1.7 % with a uniformity of the response to 
electrons of 1% and a good linearity in the energy range from 10 to 100 GeV. Th e electrom agnetic 
shower position resolution was found to be described by 1.5 mm © 5.3 mm /^/E{GeY). For an 
electron identification efficiency of 90% a hadron rejection factor of > 600 was obtained. 

PACS numbers: 29.40.Vj; 29.85.Ca; 29.85.Fj; 07.05Fb 
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INTRODUCTION 

ALICE (A Large Ion Collider Experiment) at the LHC 
is designed to carry out comprehensive measurements of 
high energy nucleus-nucleus collisions, in order to study 
QCD matter under extreme conditions and to explore 
the phase transition between confined matter and the 
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Quark-Gluon Plasma (QGP) [i|,[2]. 

ALICE contains a wide array of detector systems for 
measuring hadrons, leptons and photons. The ALICE 
detector is described in detail in 0]. The large accep- 
tance Electromagnetic Calorimeter (EMCal), which is 
foreseen to be fully installed in 2011, significantly en- 
hances Alice's capabilities for jet measurements. The 
ALICE EMCal is designed to provide the following func- 
tions: 

- efficient and unbiased fast level LO and LI trigger on 
high energy jets 

- measurement of the neutral portion of jet energy 

- improvement of jet energy resolution 



TABLE L EMCal module physical parameters. Here, RL 
stands for Radiation Length and MR for the Moliere Radius. 



measurement of high momentum photons. 



and 



electrons 

- 7 / tt" discrimination up to 30 GeV ^ 

- e / h separation (for momenta larger than 10 GeV/c) 

- high uniformity of response for isolated electromag- 
netic clusters. 

From Monte Carlo simulations, a detector energy res- 
olution of the order of about 15%/^£'(GcV) 2% was 
found to be sufficient for the jet physics program and is 
fixed as the minimum detector requirement. The elec- 
tron and photon physics programs, however, would ben- 
efit from better resolution. 

The overall design of the EMCal is heavily influenced 
by its integration within the ALICE Q setup which con- 
strains the detector acceptance to a region of about 110 
degrees in azimuth 0, —0.7 < 77 < 0.7 in pseudo-rapidity 
and 4.35 m < REMCai < 4.7 m radial distance . 

This paper presents the performance of proto- 
type modules studied in test beam measurements at 
FNAL and at CERN. The goals of these measurements 
were the determination of the intrinsic energy and 
position resolution, the investigation of the linearity and 
uniformity of the detector response, the determination 
of the light yield per unit of energy deposit and a study 
of the response to electrons and hadrons. Furthermore, 
monitoring and calibration tools were successfully 
implemented and tested. 



CALORIMETER MODULE DESIGN AND 
READOUT 

The chosen technology is a layered Lead(Pb)- 
Scintillator(Scint) sampling calorimeter with wavelength 
shifting (WLS) fibers that run longitudinally through the 
Pb/Scint stack providing light collection (Shashlik) 
The basic building block is a module consisting of 2 x 2 
optically isolated towers which are read out individually; 



Parameter 



Value 



Tower Size (at 7?=0) 
Tower Size 
Sampling Ratio 
Layers 
Scintillator 

Absorber 
Effective RL Xo 
Effective MR Rm 
Effective Density 
Sampling Fraction 
Radiation Length 



-6.0 X ~6.0 X 24.6 cm'' 
A(j)X Ar]^ 0.0143 x 0.0143 
1.44 mm Pb / 1.76 mm Scint. 
77 

Polystyrene (BASF143E + 

1.5%pTP + 0.04%POPOP) 

natural Lead 

12.3 mm 

3.20 cm 

5.68 g/cm^ 

1/10.5 

20.1 



^ Considering invariant mass and shower shape techniques only. 



each spans Arj x A(f> ~ 0.014 x 0.014 each. White, acid 
free, bond paper serves as a diffuse reflector on the scintil- 
lator surfaces and provides friction between layers. The 
scintillator edges are treated with Ti02 loaded reflector 
to improve the transverse optical uniformity within a sin- 
gle tower and to provide tower to tower optical isolation 
better than 99%. 

The energy resolution for a sampling electromagnetic 
calorimeter varies with the sampling frequency approxi- 
mately as y^dsc/fs, where dsc is the scintillator thickness 
in mm and fs is the sampling fraction for minimum ioniz- 
ing particles (MIPs). For optimum resolution in a given 
physical space and total radiation length, there is thus 
a desire to have the highest possible sampling frequency. 
Practical considerations, including the cost of the total 
assembly labour, suggest reducing the total number of 
Pb/Scint layers thus decreasing the sampling frequency. 
The requirement of a compact detector consistent with 
the EMCal integration volume and the chosen detector 
thickness of about 20 radiation lengths, results in a lead 
to scintillator ratio by volume of about 1:1.22 correspond- 
ing to a sampling geometry of Pb(1.44 mm) /Scint (1.76 
mm) . 

The physical characteristics of the EMCal modules are 
summarized in Table ID 

Scintillation photons produced in each tower are cap- 
tured by an array of 36 Kuraray Y-11 (200 M), double 
clad, wavelength shifting (WLS) fibers. Each fiber within 
a given tower terminates in an aluminized mirror at the 
front face of the module and is integrated into a pol- 
ished, circular group of 36 fibers at the photo sensor end 
at the back of the module. The 6.8 mm diameter fiber 
bundle from a given tower connects to the Avalanche 
Photodiode (APD) through a short light guide/diffuser. 
The selected photo sensor is the Hamamatsu S8664-55 
avalanche photodiode chosen for operation in the high 
field inside the ALICE magnet. The APDs are operated 
at moderate gain for low noise and high gain stability 
in order to maximize energy and timing resolution. The 
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number of primary electrons generated in the APD by an 
electromagnetic shower is ~ 4.4 electrons/McV. The re- 
verse bias voltage of the APDs are individually controlled 
to provide an electron multiplication factor (M) of 30 re- 
sulting in a charge output of « 132 electrons/MeV from 
the APDs. All APDs used for the test beam measure- 
ments were previously calibrated The charge output 
from the APD is integrated by a Charge Sensitive Pream- 
plifier (CSP) with a short rise time of « 10 ns and a long 
decay time of « 130 fxs, i.e., approximately a step pulse. 
The amplitude of the step pulse is proportional to the 
number of integrated electrons from the APD and there- 
fore proportional to the energy of the incident particle. 
The output from the CSP is conditioned with a second 
order gaussian shapcr in order to make the signal suitable 
for digitization with the Alice TPC Readout Chip 

The readout electronics of the PHOS (PHoton Spec- 
trometer) detector Q of ALICE have been adopted for 
the EMCal front end electronics readout with only minor 
modifications, as the light yield per unit of energy deposit 
in the EMCal is similar to that of the PHOS A de- 
tailed description of the EMCal (PHOS) front end elec- 
tronics (FEE) and their performance is given in Ref. 
The FEE has an effective 14-bit dynamic range over the 
interval 16 MeV to 250 GeV resulting in a Least Signifi- 
cant Bit on the low gain range of 250 MeV (10-bits) and 
on the high gain range of 16 MeV. Compared to PHOS, 
the coarse granularity of the EMCal yields higher occu- 
pancies. As the number of read out samples recorded is 
dictated by the total shaped pulse width, a shaping time 
of 200 ns (2 fis for PHOS) is chosen in order to keep the 
total data volume per central unit similar to PHOS and 
to fulfill the constraints from the total available band- 
width. This results in an electronic noise contribution of 
about 12 MeV per EMCal tower. However, due to the 
larger intrinsic energy resolution term of EMCal com- 
pared to PHOS the importance of the electronics noise 
contribution is much less. The effect of the shaping time 
on the calorimeter resolution has been studied in the test 
beam measurements performed at FNAL and is discussed 
in the next section. 



bcr 2005 at the Meson Test Beam (MTEST) at FNAL 
utilizing a stacked 4x4 array of prototype EMCal mod- 
ules (8x8 towers) of slightly different design than the 
final one, such as a radiation length of 22 Xq and a sam- 
pling geometry of Pb(1.6 mm)/Scint(1.6 mm). For this 
test in particular, measurements were made for compari- 
son of the performance with two different signal shaping 
times in the front end electronics. Two front end elec- 
tronics cards (32 towers each) were used for the readout 
of the modules; one had the nominal 2 ^s signal shaping 
time of the PHOS, and the other had a modified 200 ns 
shaping time as planned for EMCal. 

The goals of the test beam measurements were: 

• To determine the intrinsic energy and the position 
resolution using electron beams. 

• To investigate the linearity and uniformity of the 
response; in particular across towers and module 
boundaries and for tilted or recessed modules. 

• To determine the light yield (signal) per unit of 
deposited electromagnetic energy. 

• To study the effect of shorter shaping times as 
planned for the final design. 

• To study the energy dependence of the response 
to electrons and hadrons to determine the particle 
identification capabilities of the EMCal. 

• To develop and investigate the performance of the 
monitoring and calibration tools (gain stability, 
time dependencies) using electron beams, MIPs 
from hadron beams, LED events and cosmic muons. 

• To develop and test ALICE standard software for 
readout, calibration and analysis. 

Test setup and beam line instrumentation 

The characteristics of the test beams at FNAL and 
CERN are summarised in Table [III 



TEST BEAM MEASUREMENTS 

The performance of the first ALICE EMCal mod- 
ules constructed according to final design was studied 
in CERN SPS and PS test beam lines in autumn 2007. 
The test utilized a stacked 4x4 array of EMCal modules 
(8x8 towers). All towers were instrumented with the 
full electronics chain with shapers and APD gains oper- 
ated as planned in ALICE. A LED calibration system 
was installed in order to monitor time-dependent gain 
changes. The readout of the front end electronics used 
the standard ALICE data aquisition system. 

Earlier test measurements were performed in Novem- 



TABLE II: Test beam parameters. 



Lab 


FNAL 


CERN 


CERN 


test beam 


MT6 


SPS H6 


PS TIO 


particle 


e,h 


e,h 


e,h 


intensity [s~"^] 


lO^-lO"* 


10^-10^ 


10^-10^ 


Ap/p 


± 1% 


± 1.3% 




Prangc [GeV] 


3-33 


5-100 


0.5-6.5 


purity 


mixed beam 


> 99% 


mixed beam 



For handling and stacking purposes, the modules were 
assembled on a strong-back in strip units of four modules 
in the vertical direction. In order to scan the entire sur- 
face of all four modules they were placed on a remotely 
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FIG. 2: LED signal amplitude as a function of the measured 
temperature. 



were also recorded via the CCUSB system. The EM- 
Cal data were combined with the data from the trigger 
detectors and from the MWPCs offline, aligning the in- 
formation from the different data streams spill-by-spill. 



LED calibration system 



FIG. 1: Upper panel: temperature measurement as a function 
of time. Mid panel: LED signal amplitude for a typical tower 
for the same time interval. Lower panel: corrected LED signal 
amplitude. 



controlled movable platform. The range of both hori- 
zontal and vertical adjustment allowed to scan the whole 
array of modules. 

The EMCal readout electronics were attached to the 
back of the array of modules with the electronics cards 
and readout units located on the same moveable table 
as the modules, together with the low voltage supplies. 
In both setups at CERN and FNAL, a pair of scintil- 
lator paddles upstream of the EMCal was used for the 
beam definition trigger. In addition, at the CERN PS 
and at the FNAL MTEST, the signals from gas thresh- 
old Ccrcnkov counters were used as an electron trigger 
for clcctron/pion discrimination. A set of three MWPCs 
in front of the EMCal provided x — y position measure- 
ments with better than 1 mm position resolution for the 
setup at FNAL. The MWPCs were used to define the 
beam particle trajectory which could then be projected 
to the front face of the EMCal modules. 

The official ALICE data aquisition (DATE v6.13) ^ 
was used for taking the EMCal data. The MWPC data 
was recorded with a CamacCrate-via-USB (CCUSB) 
readout system. The data from the Cerenkov counters 



In order to reach the design EMCal energy resolu- 
tion for high energy electromagnetic showers, a tower-by- 
tower relative energy calibration of better than 1% has to 
be obtained and maintained in the offline analysis. In ad- 
dition, since analog tower energy sums provide the basis 
of the level LO and LI high energy shower trigger input to 
the ALICE trigger decision, the EMCal should operate 
with APD gains adjusted to match online relative tower 
energy calibrations to better than 5%. 

A LED calibration system, in which all towers view 
a calibrated pulsed LED light source, has been success- 
fully tested to track and adjust for the temperature de- 
pendence of the APD gains during operation. The LED 
triggers were collected in parallel with the beam particle 
events throughout the entire CERN test beam measure- 
ments. These measurements were performed with the 
APDs operated at the nominal fixed M=30 gain. 

The variation of the EMCal response to the LED signal 
with time and temperature was studied in order to test 
the system for calibration purposes. The temperature 
was monitored by a total of eight temperature sensors 
installed on the back surface of the module. The mea- 
sured LED signal amplitude variation for a given tower as 
a function of time is compared in Fig. [1] for the same time 
interval with the temperature readings from the nearest 
sensor for the module in which the tower was located. A 
clear anti-correlation is observed. 

Over the whole data taking period, some sharp vari- 
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ations in the LED signal amplitude were observed that 
cannot be attributed to temperature changes but rather 
to LED light yield changes, as when the setup was re- 
configured. These changes of the overall LED light were 
taken into account with an iterative extraction of the 
temperature coefficients. First, a new time interval was 
defined if an APD amplitude changed by more than 20 % 
from one hour to the next. For each time interval, both 
low and high gain LED signal amplitudes were fit simul- 
taneously as planes in space defined by signal amplitude, 
temperature, and the time interval. In a first iteration, 
all points deviating by more than 1.5 cr from a predefined 
slope range (0.015 < |dA//dT| [%/°C] < 0.025) were ex- 
cluded. In the next iteration, the cleaned sample was fit 
with a free parameter for the slope in order to define the 
temperature coefficient. Fig.[2]shows the LED amplitude 
for a typical tower as a function of the temperature and 
for a certain time interval. The temperature coefficients 
obtained from the fits of these distributions were used to 
correct for the time dependence of the APD gain. As an 
example, the corrected LED amplitude is shown in the 
lower panel of Fig. [T] for the considered time interval. 

The selected LED event amplitudes as well as the in- 
formation from the temperature sensors as a function of 
time are stored in a database. An interface was devel- 
oped and tested that allows for time-dependent calibra- 
tion corrections in the offline analysis of the test beam 
data. 

Signal reconstruction 

The digitized time samples from the read out have an 
amplitude as a function of time t that can be described 
with the form of a F-fimction in ADC{t), where 

ADC{t) = pedestal A • e"" • • e" (i"^) , (1) 

X = {t- to)/T . 

Here, r = n • tq with the shaper constant tq and n=2 as 
the shaper is gaussian of second order (composed of a dif- 
ferentiator and two integrators The charge collected 
from the APD, and hence the energy deposited in the 
tower, is proportional to the value of the parameter A at 
the time value {to + t) where the function peaks. 

The test beam data were used to investigate the per- 
formance of this function and the parameters were opti- 
mized. The High-Low gain correlation was studied using 
the electron data in order to determine a threshold value 
for the amplitude for which the low gain rather than the 
high gain needs to be used due to saturation (at 1023 
ADC counts). A good High-Low gain correlation with 
an average ratio of 16.3 between both gains was found 
up to 1050 ADC counts when using the values from a fit 
for ADC counts > 1000. 

An overall inter-calibration procedure was carried out 
for all towers by normalizing the hadron MIP amplitudes 
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FIG. 3: Reconstructed energy for 80 GeV incident electrons. 
The curve represents a fit of a truncated gaussian to the his- 
togram with fit results as given in the figure. 

in each tower, to a reference tower. Isolation of the MIP 
peak was achieved requiring, for each tower, no energy 
deposit in the surrounding eight towers. An alternative 
inter-calibration map was also considered by using the in- 
formation given by the electron beam peak in each tower. 
An absolute calibration for each tower was accomplished 
by comparing the nominal electron beam energy with the 
corresponding peak in the energy spectrum, as obtained 
by a sum over a 3 x 3 tower cluster. For this purpose, 3x3 
local cluster inter-calibration coefficients were extracted 
from the overall map, by choosing each tower in turn as 
a reference. This allowed evaluation of the energy spec- 
trum by a sum over the 9 towers in the cluster, with a 
proper calibration adjusted to match that of the central 
tower in each cluster. These calibration coefficients were 
used to analyse the test beam data with the standard 
ALICE cluster reconstruction software. Fig. [3] shows the 
reconstructed energy for 80 GeV incident electrons (for a 
typical run). 

Linearity and uniformity of energy response 

The absolute energy calibration obtained by a sum 
over a 3x3 tower cluster is shown in Fig. (4] The lin- 
earity of the response is better than 1% over the full 
energy range down to 20 GeV. At low energies, thresh- 
old effects become non-negligible compared to the to- 
tal energy deposited and light transmission losses might 
have an impact. In fact, as shown in the bottom panel 
of Fig. 3] by the full circles the reconstructed energy is 
systematically lower than the incident one for energies 
equal or below 10 GeV. A drop of ~ 10% is observed 
at 5 GeV. This behaviour is well described by a cubic 
function as demonstrated by the open circles in Fig. [4] 
(bottom panel). At high energies, deviations of the ratio 
from unity are expected due to longitudinal shower leak- 
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it is also important to estimate the tower-to-tower dis- 
persion of the Hght yield. 

During the test beam, the APDs were all operating at 
gain M=30. The individual voltage settings had been 
established for each APD prior to the test beam mea- 
surements. This procedure compares the amplitude at 
a given bias voltage to the amplitude measured at low 
voltages, where the gain is assumed to be unity The 
light yield (LY) at the gain M=30, for each individual 
tower, is then extracted following 

LY{p.e./MeV) = {channels / M eV) ■ {1/G a) (2) 

•(1/Pg) • {l/ADC,onv), 

where the shaper amplifier gain Ga = 0.229, the charge 
voltage conversion factor of the preamplifier Pq = 0.83 
V/pC and the ADC conversion ADCconv — 1024 chan- 
nels/ V. The light yield at unit gain (M=l) is obtained 
from this value divided by 30. An average light yield of 
about (4.3 ± 0.3) photoelectrons/MeV, was found, which 
is consistent with the light yield value of PHOS. 

Energy resolution 

The energy resolution of an electromagnetic calorime- 
ter can be parameterized as 



FIG. 4: Top: Linearity of the response for a sum over a 3x3 
tower cluster as a function of the incident beam energy. The 
solid line is a linear fit to the data. Bottom: Ratios between 
the linear fit and the data (full circles) and a cubic fit and 
the data (open circles). The dashed line is placed at a ratio 
of unity to guide the eye. 

age. The data show an indication of such an energy loss 
at high energies. 

The uniformity of the energy response was studied for 
several different conditions. All module centers and a 
major part of tower centers were scanned using 80 GeV 
electrons. In addition, data were taken across tower and 
module borders. A uniformity of the energy response 
was found with a RMS better than 1 GeV, for 80 GeV 
incoming electrons. This result implies a very good uni- 
formity of the response (within 1%) for the EMCal as 
constructed. 



Light Yield 

The light yield, the number of photoelectrons at the 
APD per unit of electromagnetic energy deposited in 
the EMCal (photoelectrons/MeV), determines the over- 
all APD-l-shaper gain required to match the desired dy- 
namic range in ALICE. Due to the large number of indi- 
vidual towers planned for the final design of the EMCal, 



(t{E)/E = a®b/y/E®c/E , (3) 

where E is the measured energy. The intrinsic resolu- 
tion is characterized by the parameter b that arises from 
stochastic fluctuations due to intrinsic detector effects 
such as energy deposit, energy sampling, light collection 
efficiency, etc. The constant term, a, originates from 
systematic effects, such as shower leakage, detector non- 
uniformity or channel- by-channel calibration errors. The 
third term, c, is due to electronic noise summed over 
the towers of the cluster used to reconstruct the electro- 
magnetic shower. The three resolution contributions add 
together in quadrature. 

Detailed GEANT3 Monte Carlo simulations for the 
final module design yield fit results using Eqn.Q of 
a = (1.65±0.04)%, b = (8.0±0.2)% and c = (7.4±0.2)% 
over a photon energy range of 0.5 GeV to 200 GeV. 
These results are based on energy deposition only and 
do not include photon transport efficiencies. Systematic 
contributions to the resolution arising from calibration 
and related systematic uncertainties are ignored. The 
value of the constant term a is dominated by longitudi- 
nal shower leakage in these calculations. Other system- 
atic effects, which arise during detector fabrication and 
from the tower- by-tower calibration uncertainties, will in- 
crease a. 

By combining data taken at the CERN PS and SPS the 
calorimeter resolution over the energy range of 0.5 GeV 
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FIG. 5: Energy resolution for electrons as a function of the 
incident beam momentum. The beam energy spread was sub- 
tracted from the measured result. The dashed curve repre- 
sents the resolution obtained from Monte Carlo simulations. 
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FIG. 6: Energy resolution for electrons as a function of the in- 
cident beam momentum for short (full circles) and long (open 
circles) shaper time, corresponding to the EMCal and PHOS 
design, respectively. The open circles are slightly shifted to 
the right for visibility. 



to 100 GcV could be explored. Such energy scans were 
performed at several different positions, including tower 
and module edges. The LED calibration system was used 
to track and adjust for the time dependence of the cal- 
ibration coefficients. No systematic variation of the res- 
olution depending on the position was observed. The 
resolution obtained at the different positions was com- 
bined and the average values as a function of the incident 
beam momentum arc displayed in Fig. [5l The momen- 
tum spread of the incident beam of typically 1.3% was 
subtracted in quadrature. A fit to the energy resolution 
as a function of the incident energy following Eqn.([3]) is 
also shown in Fig. [5] with parameters [a = 1.7 ± 0.3)%, 
(b = ll.l±0.4)%and (c = 5.1±0.3)%. These parameters 
can be compared with the GEANT3 simulation result for 
the EMCal module geometry described before and pre- 
sented by the dashed line in Fig. O The increase of the 
stochastic term 6, representing a worse intrinsic resolu- 
tion compared to the Monte Carlo simulations, is mainly 
due to light attenuation and light collection inefficiencies 
which were not modelled. The small increase of the con- 
stant term a demonstrates a stable, high quality detector 
fabrication and a good tower-by-tower calibration. The 
linear term, modelling electronic noise contributions, is 
set too high in the simulation. 

The energy resolution was also studied for different 
incidence locations corresponding to the modules as in- 
stalled in ALICE. Most of the test beam data were taken 
with a configuration where the beam hits the EMCal 
modules perpendicularly, corresponding to z = 0, 77=0 
position. Data were also taken with configurations 
where the modules were tilted in by 6 or 9 degrees 
at different surface positions. The energy resolution for 
such tilted configurations compares well with the average 



resolution as a function of energy presented in Fig. [S] 
No significant deviations from the average resolution at 
zero degree was observed. 

Using the data from the FNAL test beam, possible 
effects of the shorter design shaping time for the EMCal 
of 200 ns (compared to 2 fis for PHOS) were studied. 
Fig. [6] shows the energy resolution as a function of the 
incident energy. The results are shown separately for the 
short (full circles) and long (open circles) shaping time 
readout regions of the test setup, averaged over various 
runs in each region. The resolution slightly deteriorates 
when using the short shaping time but is still well within 
the detector requirements. 



Position resolution 

The segmentation of the calorimeter allows one to ob- 
tain the hit position from the energy distribution inside 
a cluster with an accuracy better than the tower size. 
The X and y coordinate locations are calculated using a 
logarithmic weighting ll| of the tower energy deposits. 
Data from the FNAL test beam were used where the 
MWPCs provided a reference position measurement of 
better than 1mm. Fig. [7] shows the x and y position 
resolution as a function of the energy deposit for elec- 
trons. As expected, no difference in the resolution in 
the X and y directions is observed. The electromagnetic 
shower position resolution is seen to be described as 1.5 
mm © 5.3 mm/ ^y E(GeV) , where the two contributions 



add together in quadrature. 
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FIG. 7: Dependence of the position resolution as a function 
of the deposited energy for electrons. The curve shows the 
best fit result. The triangles representing the resolution for 
the y-position are slightly shifted to the right for visibility. 



Response to hadrons 

The EMCal can further enhance the ALICE parti- 
cle identification capabihties due to the characteristically 
different response to electrons and hadrons. While elec- 
trons leave all their energy in the calorimeter, hadrons 
leave only a fraction of their energy but show a long tail 
due to hadronic showers. At the CERN PS pure elec- 
tron and hadron beams were available. Fig. [8] shows the 
reconstructed energy for an electron and hadron beam 
of 100 GeV, which illustrates this very distinct response 
to electrons and hadrons. The high energy tail in the 
hadron response originates from processes such as charge 
exchange tt^ + p —^ tt^ + n, where most of the energy 
of the charged pions goes into neutral pions. These neu- 
tral pions decay immediately into photons starting a cas- 
cade which is indistinguishable from an electron-initiated 
shower. 

The hadron rejection factor is defined as the number of 
all hadrons divided by the number of hadrons misiden- 
tificd as electrons. This factor is shown in Fig. [9] as a 
function of the incident hadron beam energy for electron 
identification efficiencies of 90% and 95%. Error bars give 
the total uncertainty, which is dominated by the system- 
atic uncertainty of the evaluation. Results from a Monte 
Carlo simulation arc also shown for an electron identifica- 
tion efficiencies of 90%. For each incident beam energy, 
the electron identification efficiency was determined by 
integrating the reconstructed energy distribution of the 
pure electron beam (dashed histogram in Fig. [5]) from 
the right-hand side till a cut value corresponding to the 
chosen efficiency. 

A rejection factor of 10^ to 10^ is obtained over the 
energy range of 40 GeV to 100 GeV. Test beam data at 
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FIG. 8: EMCal response to hadrons (full histogram) and elec- 
trons (dashed histogram) of 100 GeV. 
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FIG. 9: Hadron rejection as a function of the incident hadron 
beam energy for an electron identification efficiency of 90% 
(circles) and 95% (squares). Error bars represent the to- 
tal uncertainty. The open circles show the result from a 
Monte Carlo simulation for 90% electron identification effi- 
ciency. The squares (open circles) are shifted to the left (right) 
for visibility. 



lower hadron energies were not taken. Hadron/electron 
rejection can be further improved by considering the 
characteristic shower shapes, as hadrons produce showers 
with wider spatial distributions than electrons. 



Cosmic ray calibration 

A calibration of all modules will be performed before 
their insertion in ALICE. The calibration procedure is 
based on a measurement of cosmic-ray muons at the min- 
imum of ionization. 

The muon signal measured in each tower is obtained by 
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FIG. 10: Response of 384 towers of the EMCal to cosmic-ray 
muons before (dashed histogram) and after (full histogram) 
individual gain calibration. The curve represents a fit of a 
gaussian to the full histogram with fit results as given in the 
figure. 



the use of an isolation procedure applied offline. For each 
event, the maximum signal amplitude is chosen and for 
all neighboring towers a signal smaller than a threshold 
value required. This threshold value is limited by the 
electronic noise (set to 3 ADC channels in the present 
case which amounts to about 15% of the muon energy). 

Since the energy of MIP muons is too low to trigger 
the EMCal, an external trigger is necessary. The muons, 
passing the towers along their length, are selected us- 
ing scintillator paddles. Each paddle covers 12 modules 
grouped into a 'strip module', and is read out at both 
extremities by photomultiplier tubes (PMT). This trig- 
ger configuration appeared to be the most reliable from 
the cosmic analysis test done in December 2007 Q with 
the EMCal prototype described above. The time of flight 
difference between both PMTs allows one to select verti- 
cal muons with a spatial accuracy of a few centimeters. 
The isolation procedure then ensures that no energy was 
deposited in the neighboring towers. A 24- hour run al- 
lows the accumulation of about 500 muons per towers, 
which is sufficient to extract a MIP peak with an accu- 
racy better than 1%. 

An individual gain calibration is performed for each 
tower, so as to ensure that the amplitude of the aver- 
age signal for cosmic muons is the same for all towers. 
The tower gains, which are controlled through the tower 
high voltage power, are tuned iteratively. Fig. [10] shows 
the dispersion of the mean amplitude of 384 towers be- 
fore and after this procedure (thin and bold lines, respec- 
tively). After three iterations a final relative dispersion 
< 3% is reached. 



CONCLUSION 

The performance of a (4 x 4) array of prototype mod- 
ules and of a (4 X 4) array of final design modules for the 
ALICE EMCal has been studied in test beam measure- 
ments at FNAL and CERN, respectively. 
These studies demonstrate: (i) an average light yield of 
(4.3 ± 0.3) photoelectrons/MeV 

(ii) an energy resolution of 



(1.7±0.3)' 



(11.1 ±0.4) (5.1 ±0.7) 

VS(GeV) ® E{GeV) 



(iii) a uniformity of the response within 1% for all tow- 
ers and configurations 

(iv) a good linearity of the response to electrons in the 
energy range 10-100 GeV 

(v) an only slightly deteriorated energy resolution 
when using the EMCal default shaping time of 200 ns 
compared to 2 for PHOS 

(vi) a position resolution described by 1.5 mm © 5.3 
mm/^£;(GeV) 

(vii) a hadron rejection factor > 600 for an electron 
identification efhciency of 90%. 

A LED calibration system was successfully tested 
to track and adjust for temperature dependent effects 
during operation. 

Cosmic ray calibrations allow a precalibration of all 
modules prior to installation in ALICE with a relative 
spread of < 3%, sufficent for the use in an online trigger. 

The performance of the tested EMCal modules reaches 
all design criteria. 
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